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Abstract 

In this article we present a high-precision evaluation of the expansions in e = (4 —d)/2 of (up to) four- 
loop scalar vacuum master integrals, using the method of difference equations developed by S. Laporta. 
We cover the complete set of ‘QED-type’ master integrals, i.e. those with a single mass scale only (i.e. 
rrii £ {0, m}) and an even number of massive lines at each vertex. Furthermore, we collect all that is 
known analytically about four-loop ‘QED-type’ masters, as well as about all single-mass-scale vacuum 
integrals at one-, two- and three-loop order. 


1 Introduction 

Higher-order perturbative computations have become a necessity in many areas of theoretical physics, be 
it for high-precision tests of QED, QCD and the standard model, or for studying critical phenomena in 
condensed matter systems. 

Most recent investigations employ a highly automated approach, utilizing algorithms that can be im¬ 
plemented on computer algebra systems, in order to handle the growing numbers of diagrams as well as 
integrals which occur at higher loop orders. 

Computations can be divided into four key steps. First, the complete set of diagrams including symmetry 
factors has to be generated. For a detailed description of an algorithm for this step for the case of vacuum 
topologies, see Ref. pp. Second, after specifying the Feynman rules, the color- and Lorentz-algebra has 
to be worked out. Third, within dimensional regularization, massive use of the integration-by-parts (IBP) 
technique [2] to derive linear relations between different Feynman integrals in conjunction with an ordering 
prescription [3] can be used to reduce the (typically large number of) integrals to a basis of (typically a 
few) master integrals. Practical notes as well as a classification of vacuum master integrals are given in 
Ref. 0 . Fourth, the master integrals have to be solved, either fully analytically, or in an expansion around 
the space-time dimension d of interest. It is the fourth step that we wish to address here. 

A very important subset of master integrals are fully massive vacuum (bubble) integrals, since they 
constitute a main building block in asymptotic expansions (see e.g. Ref. [5]). They are also useful for 
massless theories, when a propagator mass is introduced as an intermediate infrared regulator |^. In four 
dimensions, this class of master integrals has been given up to the 4-loop level in Ref. [7]. As an application, 
these integrals are vital for computing the 4-loop QCD beta-function and anomalous dimensions |S]. In 
lower dimensions, perturbative results are needed for applications in condensed matter systems, as well as 
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in the framework of dimensionally reduced effective field theories for thermal QCD, where recent efforts 
have made four-loop contributions an issue [H] ■ We have recently extended the work of Ref. [7] , to give the 
complete set of fully massive vacuum master integrals in three dimensions, again up to the 4-loop level cni. 

The next larger set of scalar vacuum master integrals are those in which there is only one mass-scale to, 
i.e. the propagators l/(pf -I- to^) have masses rm G {0,to}. These integrals are needed for problems with 
widely separated mass scales, in which one then sets the masses of all heavy particles to to and those of all 
light particles to zero. As a well-defined subset of these single-mass-scale integrals, we here treat ‘QED-type’ 
vacuum integrals, i.e. those with an even number of massive lines at each vertex, at the 4-loop level. A 
recent application is in the computation of heavy-quark vacuum polarization m 

The complete set of ‘QED-type’ vacuum master integrals up to the 4-loop level has already been identified 
in Ref. ^ . The main purpose of this work is to numerically compute this set in terms of a high-precision e- 
expansion in d = 4—2e dimensions, and to present new analytic results for some low-order (in e) coefficients. 
Eurthermore, we have made an attempt to collect all presently known analytic results on 4d single-mass-scale 
vacuum integrals, up to four loops, in a coherent notation. 

The plan of the paper is as follows. In Sectional we give a brief review of the method of difference 
equations applied to vacuum integrals. In Section|Sl we discuss the actual implementation of the algorithm. 
In Sectional we display our numerical results for the truncated power series expansions in e of our master 
integrals, up to the four-loop level, in d = 4 — 2e. In Section [3 we discuss one case of a master integral 
which we needed to solve via a Laplace transform of its difference equation. In Section we list analytic 
results. 


2 The evaluation of master integrals through difference equations 

The method we have chosen to compute the coefficients of the truncated power series expansions of the 
master integrals is based on constructing difference equations for the integrals and then solving them nu¬ 
merically using factorial series. This approach has recently been developed in Ref. ^j, and below we briefly 
summarize its basic concepts following the notation of the original paper, which contains a much more 
detailed presentation of the subject. While the method is completely general as it applies to arbitrary kine¬ 
matics, masses and topologies na, our brief summary is somewhat adapted to the specific case of vacuum 
integrals. 

The main idea is to attach an arbitrary power x to one of the massive^ lines of a master integral U, 

" I Dfih^ ■ 

where the Di = {pf -|- ml) denote inverse scalar propagators. In our case the mass parameter mi has only 
two values, 0 and to, the latter of which we set to 1, noting that it can be restored in the end as a trivial 
dimensional pre-factor of each integral. The original master integral is then just U = D(l). Depending on 
the symmetry properties of the integral, there can be different choices for the ‘special’ line with the arbitrary 
power X, but in the limit a: = 1 they all reduce to the original master integral U. This degeneracy can (and 
will later) be used for non-trivial checks of the method. 

Employing IBP identities in a systematic way, it is possible to derive a linear difference equation obeyed 
by the generalized master integral U{x), 

R 

'^Pj{x)U{x + j) = F{x) , (2.2) 

_ 

^The massiveness is crucial in order to avoid problems with the infrared behavior of the integral. 
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where i? is a finite positive integer and the coefficients pj are polynomials in x (and the space-time dimension 
d). The function F on the r.h.s. is a linear combination of functions analogous to U{x) but derived from 
‘simpler’ master integrals, i.e. integrals containing a smaller number of loops and/or propagators. 

The general solution of this kind of an equation is the sum of a special solution to the full equation, Uq{x), 
and the solutions to the homogeneous equation (F = 0), 


U{x) 


Uq{x) + ^Uj{x), 


f=i 


where each {j = 0,i?) 


U,{x) = F 


X \ ^ 

'i ' 

s=0 




r{x + i) 


r(a; -I- 1 -I- s — Kj) 


(2.3) 


(2.4) 


is a factorial series^. Substituting this form into Eq. lOJ, one obtains the coefficients pj and Kj (the latter 
being a function of d ), as well as recursion relations for the x-independent coefficients aj (s) (being functions 
of d as well) for each solution. For the homogeneous solutions, these recursion relations relate all coefficients 
with s > 0 to their (in principle arbitrary) value at s = 0, aj{s) = Cj(s) aj(0), where the Cj{s) are rational 
functions of d. For the special solution, all ao(s) are on the other hand completely fixed in terms of the 
inhomogeneous part F'(x), consisting of ‘simpler’ integrals which are assumed to be already known in terms 
of their factorial series expansions. 

What clearly remains to be done is to fix the x- and s-independent constants 0^(0), j ^ 0, in order to 
determine the weights of the different homogeneous solutions. To this end, it is most useful to study the 
behavior of U{x) at large x. Writing the integral in the form 

= / JpfTW 

it is easy to see that its large-x behavior is determined by the small-momentum expansion of the two-point 
function which has one loop less than the original vacuum integral. 

In the case of integrals for which the limit g(0) is well-defined and non-zero, the calculation becomes 
particularly simple. Then the large-x limit of U(x) factorizes into a one-loop bubble carrying the large 
power X and a lower-loop vacuum bubble g(0), which corresponds to I7(x) with its ‘special’ line cut away. 


lim U(x) 

X —»-oo 



1 

{pI + 1)" . 


X 



{ITx-'^/^glS)) . 


( 2 . 6 ) 


A comparison with the large-x behavior of Eqs. JOll . itTHl . proportional to Ylj : can now be 

used to fix the aj(0), of which maximally one will turn out to be non-zero for our set of integrals. 

If on the other hand g(0) = 0, the treatment of the small-pi limit of this function becomes more involved. 
Fortunately, the massless lines of the sub-diagram — which were responsible for the vanishing of its value 
at zero external momentum in the first place — also make its analytic evaluation more straightforward. 
Performing a careful analysis of the subgraph, one always ends up with an integral of the type 


lim C/(x) 

X—^CXD 


/ 


(p^r 

(pf +1)’” 


(2.7) 


from which the calculation proceeds just as above providing us with the values of the 0 ^( 0 ). 

Having the full solution at hand, we have in principle completed our task, as in the limit x = 1 we recover 
from C/(x) the value of the initial integral. Let us, however, add a couple of practical remarks here. What 

■^For a rigorous definition of the concept as well as a motivation for this kind of an Ansatz, we refer the reader to Ref. 1^. 
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is still to be done is to perform the summation of the factorial series of Eq. TIM . which in practice means 
truncating the infinite sum at some large but finite Smax- Studying the convergence behavior of these sums, 
one notices that even in the cases where they do converge down to x ^ 1, their convergence properties 
usually strongly decline with decreasing x. This means that in practical computations, where one aims at 
obtaining a maximal number of correct digits for 17(1) with as little CPU time as possible, the optimal 
strategy is to evaluate the integrals U{xniax + 1), •••, U(xmax + R) with the factorial series approach at some 
a^max ^ 1 and then use the recurrence relation of Eq. 12.211 to obtain the desired result at a; = 1. The 
price to pay is, however, a loss of numerical accuracy at each ‘pushdown’ {x ^ x — 1) step due to possible 
cancellations, which makes the use of a very high x^ax impossible. In practice the strategy is to determine 
an optimal value for the ratio Smax/a^max- To give an example, for the four-loop integrals of Section 0] we 
have found that Smax/a^max 20...40 is a good value, while we used a range of Smax 2000 ... 2500. Eor a 
few special cases, for which additional numerical problems emerged, we were forced to limit the value of the 
parameters to roughly Smax 200 and Xmax ^ 30, which decreased the accuracy of the results significantly. 

3 Implementation of the algorithm 

As is apparent from the preceding section, there are three main steps involved in obtaining the desired 
numerical coefficients in the e-expansion of each master integral: deriving the difference equations obeyed 
by each integral, solving them in terms of factorial series, and finally performing the e-expansion and 
numerically evaluating the sum of Eq. TIM (t runcated at Smax) to the precision needed. We will briefly 
address each of them in the following. 

Eor the first step, we slightly generalized the IBP algorithm we had used for reducing generic 4-loop bubble 
integrals to master integrals, which follows the setup given in Ref. [3], and whose implementation in EORM 
m is documented in Ref. 0. The main difference is an enlarged representation for the integrals, keeping 
track of the line which carries the extra powers x, as well as the fact that there are now two independent 
variables (d, x), requiring factorization (and inversion) of bivariate polynomials, as opposed to univariate 
polynomials in the original version. 

Second, staying within EORM for convenience, we implemented routines that straightforwardly solve 
the difference equations in terms of factorial series, along the lines of Ref. |3]- This is done starting with 
the simplest one-loop master integral, and working the way up to the most complicated (most lines) four- 
loop integral, ensuring that at each step, the ‘simpler’ terms constituting the inhomogeneous parts of the 
difference equation are already known. The output are then plain ascii files specifying each solution in the 
form of Eq. (E3I) as well as containing recursion relations for the coefficients a{s). Note that these first two 
steps are performed exactly, in d dimensions. 

Third, once the recursion relations for the coefficients a{s) were known, we used a Mathematica program 
to obtain their numerical values at each s to a predefined precision, and to perform the summation of 
the factorial series. While this procedure is in principle straightforward, there are some twists that we 
employed to help reduce the running times significantly, most of which are probably quite specific to our 
use of Mathematica. To avoid a rapid loss of significant digits in solving the recursion steps that relate each 
a{s) to a(0), especially those for the homogeneous coefficients, we first solved the relations analytically and 
only in the end substituted the numerical value (actually the truncated e-expansion) of the first non-zero 
coefficient. In fact, we found Mathematica to operate quite efficiently with operations like multiplication of 
two truncated power series, so that we relied heavily on it. Eurthermore, since — not surprisingly — the 
most time-consuming part in the summation of the series turned out to be the e-expansion of T-functions, 
we achieved a notable speed-up by substituting the T-functions with large arguments by suitable products 
of linear factors times T-functions of smaller arguments. Finally, a vital step in avoiding an excessive loss in 
the depth of the e-expansions when going from one integral to the next, was to apply the ‘Chop’ command 
to remove from the results and coefficients excess unphysical poles, whose coefficients were of the order of, 
say, 10“®'' or less. In some cases we were in addition able to reduce the loss of precision in the pushdown 
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steps by first analytically solving 17(1) as a function of 17(xniax)) and only in the very end substituting the 
numerical value of the latter. 


4 Numerical results 


Below we list the Laurent expansions in e = (4 — (i)/2 of vacuum master integrals up to four loops. We use 
an intuitive graphical notation, in which each solid line represents a massive scalar propagator l/(p^ + 1) 
and a dashed line a massless one 1/p^. The integral measure we have chosen here is 


1 


j4-2£^ 


r(2 + e) J tv' 


2-e 


(4.1) 


which implies that the 1-loop tadpole is J = — ~ ^provide 

the results to order keeping the accuracy at 50 significant digits for the 2- and 3-loop master integrals 
and at 40 for the 4-loop ones. There are two exceptions. For one of the 3-loop integrals (see Eq. (lOill 
below) the factorial series does not converge and hence the integral has to be treated by Laplace transform, 
see Sectional For one of the 4-loop integrals (see Eq. (14.1411 below) we only give the first seven e-orders to 
17 significant digits. To obtain more e-orders and significant digits for all integrals listed here is merely a 
matter of additional CPU time. 

We have produced numerical results for all single-mass-scale vacuum master integrals up to three loops 
(these are the master integrals entering the package of Avdeev jTl] and MATAD 53 )i and for all ‘QED- 
type’ vacuum master integrals at four loops. Here, we display only those numerical results which correspond 
neither to analytically solvable integrals (1 of 1 1-loop master, 1 of 2 2-loop masters, 3 of 12 3-loop masters, 
2 of 10 4-loop masters, all of which are given analytically in Section M below, and are listed in numerical 
form in the appendix), nor to fully massive cases (1 of 2 2-loop masters, 3 of 12 3-loop masters, which are 
given in Section 9.3.1 of Ref. 0; 1 of 10 4-loop masters, which is given in Eq. (4) of Ref. 5]). 



-f 1.0000000000000000000000000000000000000000000000000 

-f 0.75000000000000000000000000000000000000000000000000 e-2 
-f 2.8750000000000000000000000000000000000000000000000 
-f 1.8362912870512825038535151499626234054646567716807 

- 26.427828097688527267319254765120590367456377175480 e 

- 35.088051385481306364065961402117419432373775682177 

- 512.75537623727044689027104289971864365971796649684 

- 607.61494953927726782115473332930225595551912885034 

- 5868.5987295458313170081280447224279031237930577453 e® 

- 6835.6108788455114123641492279253803965001408075543 e® 

- 58194.090725773231428299235587057139067942816045554 

- 67435.335245041201792055506063164867635607825896649 e® 

- 546094.78026628005592146280450252032502449454982782 e® 

- 631563.41278231233491152773645513360004834876043263 + O (e^^) (4.2) 

- 0.66666666666666666666666666666666666666666666666667 e“® 

- 1.6666666666666666666666666666666666666666666666667 
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- 5.1290732088140381005700717333376095855758453867530 

- 26.359970069205366659319388577532454678949629074714 

- 27.711175418518951962132692178901111387631205211816 e 

- 293.15661097603756640443665077615751632698451158842 

- 142.70296384808301760189570443963964069968530061393 

- 2882.1838924952422595902727649575335612132315437366 e'‘ 

- 801.64629651874722343778241421866459175486305997074 

- 26947.975116190322227046885628024191588708203470044 e® 

- 5202.1954102253813787831194867097139379256134207497 

- 246612.58893683893330836716807041349553641918919392 e® 

- 38662.312198716830334636275721442722552625311137805 e® 

- 2235893.9169450155346378997842831790622571843918826 + O (e^^) (4.3) 

(X^. = + 2.4041138063191885707994763230228999815299725846810 e-i 

- 13.125546202841586242894146861604104971473328745577 
+ 58.026260003878655576719786597271170572487856789112 e 

- 215.15799420773251496754795359758001229751012774168 
+ 741.02167568175382570477503744319405056068752490840 e® 

- 2422.8745603243623464433277838674972111388328822910 
+ 7691.0946660371679072096695375419004253731139722117 e® 

- 23935.477541938694107878632636617038283240279946231e® 

+ 73567.948130076321368433008921329180208448045889323 

- 224259.22429731742234354745849077250951082319381698 e® 

+ 679949.06185664528482517935972319009053625673803897 e® 

- 2054250.6137900838709156880585181458050495843111626 + O (e^^) (4.4) 

(X) = + 2.4041138063191885707994763230228999815299725846810 

- 10.073203643096893062671213536841941862151359216063 
+ 46.082030897278984204342981632818973100797752268016 e 

- 162.84321571472549604685427998929495247564452607777 
+ 563.02541599052549921690912303391142482056503193963 e® 

- 1822.8278416039379661792322993085062379900421439244 e'‘ 

+ 5785.9815122286472118701238800303861843636492370028 e® 

- 17968.847688521304415691142884872355494614421789034 e® 

+ 55216.376506037509642111329809657803343085975716148 

- 168240.56307714987438328576061703042187348961355271 e® 

+ 510052.27830760883492002666904963035268124951366379 e® 

- 1540802.7858406456522499592944279183737583148997338 e®° + O (e®®) (4.5) 

= + 2.40411380631919 e"®- 6.09209302191832 



6 


+ 35.8130598712514 e - 104.744695525740 
+ 394.7643404810 - 1200.978166746 + O (e^) 


(4.6) 



+ 2.4041138063191885707994763230228999815299725846810 

- 10.239350912945217732184803670827657230740659540460 
+ 46.310233388509835938575195677581891346572247104081e 

- 163.71903666846274587940160817767510251267798801563 
+ 564.30910069499791449917891053192483169414448523830 

- 1825.8206691490586101339592917414250683553208417056 e'‘ 

+ 5790.4830503256226500389801844087331481909799175043 e® 

- 17977.329926014373828927297954140399343188401436964 e® 

+ 55229.228709840982549894271961274899889724333654202 

- 168262.33984881039469900336415583383277506556682677 e® 

+ 510085.27212468040781599787353617454516173782111015 e® 

- 1540855.6379207615212777547545938657889645397412471 + O (e^^) (4.7) 



1.000000000000000000000000000000000000000 


- 0.500000000000000000000000000000000000000 e"® 

- 3.527777777777777777777777777777777777778 

- 1.995370370370370370370370370370370370370 

- 36.82021604938271604938271604938271604938 

- 19.87920801451107035069575635156380101575 e 

- 1809.001126638637160933894507798781706682 

- 941.2486498215135407529753614624254521594 e® 

- 49114.80404240275263940837370626747663512 

- 25712.87944658606239888301931387195377680 e® 

- 1014742.540337323108931396794699794706304 e® 

- 533925.9315185165221824312193117157135164 

- 18513953.44519328478360685998151320048728 e® 

- 9769845.146715270007449428486953016122496 e® 

- 317669932.9515691976277658362596784695115 + O (e“) (4.8) 



+ 0.2500000000000000000000000000000000000000 
+ 0.5000000000000000000000000000000000000000 e"® 
+ 1.000000000000000000000000000000000000000 £-2 
+ 1.813369870537362855098298049824424939972 e"® 

- 113.8224542836131461311762552843948945680 

- 33.70692008121875082746730709549318292582 e 

- 3800.131177952398833364468086486701310324 

- 724.2483980459868435529785916706580415218 e® 
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- 83243.75114211351557600351242603548310943 

- 9962.244874731471054690629554209449745080 

- 1556494.392681571176934758495668112116219 e® 

- 125852.9269007094630774780949002157883896 

- 27026768.74324139691004925806420463865625 e® 

- 1619900.985945231199760429618558131115494 

- 451968203.1707264233126870326577507342793 + O (e“) 

+ 0.6666666666666666666666666666666666666667 
+ 1.333333333333333333333333333333333333333 e"® 

+ 3.333333333333333333333333333333333333333 

- 2.922363183148830477868063138605600049253 e"® 

- 52.50529739842769756973487794955803028226 

- 622.1548972708590376012515685880304077291e 

- 1741.392944346262052260405956917114927201 

- 17196.12685330902582768340098554237636824 e® 

- 35037.76438140725371856293904191777384497 

- 350040.6052285494016783912074340410365119 e® 

- 619669.7756160060500505884016704452111642 e® 

- 6316632.078794015469602341973684043315269 

- 10420684.66626276045383010214928284087492 e® 

- 107682720.9936656086807201002498313447872 e® 

- 171175743.7334785889316026497774359587050 e^® + O (e“) 

- 0.1666666666666666666666666666666666666667 

- 0.8333333333333333333333333333333333333333 e"® 

- 5.535390236492927618733071494844783324098 e"^ 

- 18.82211358179364443034084677047078519365 e"® 

- 25.33131709639103630297934297632219102642 

- 692.6253681383859207802291611352811358818 e 
+ 1304.406827189023173835521731683389467596 

- 17597.62761742767175796342110253040416842 e® 

+ 43608.68478725040973761321535022863250602 

- 356925.7947952212585233385307804100264907 e® 

+ 939175.5208936133499000732308171881559393 e® 

- 6467516.567931160982387324881460909595434 
+ 17364082.00316543469946942544036134544207 e® 

- 110630064.0504718962799294108969364410321 e® 

+ 299555848.2967199841801845664112146544429 e^° + O (e“) 
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(4.9) 


(4.10) 


(4.11) 




- 0.1666666666666666666666666666666666666667 

- 0.8333333333333333333333333333333333333333 

- 4.934361784913130476033202414089058328716 e"^ 

- 16.20728580177089935860423675579773428042 

- 66.26045267088253719950145322703822192330 

- 375.8131807648258568590885739785987842170 e 

- 558.8684980291056611497739327005648555056 

- 8005.909258052308131324607228836046381631 

- 1919.455472313357401527573620257269458394 

- 151678.4473872174906037102434226085312395 e® 

+ 45367.13703676553118943884640644928312337 e® 

- 2666691.771475554201053771443874777744895 

+ 1469487.259313641787542428480707831779115 e® 

- 44948193.38776277686896804573212402472750 e® 

+ 30458142.54983328875970186426762508169587 + O (e“) 

+ 5.184638775716849631656827432285170840285 

- 43.27615856932464061120186605318440014787 
+ 281.6878028207571441503294898731325032572 e 

- 1513.439498357334510783856768928671867962 

+ 7425.188218180801392788420314778806091908 e® 

- 34157.28328369996440389044714608216014037 
+ 150927.1863992836861076852919478736657397 e® 

- 648212.2200729992344766408701706025837197 e® 

+ 2730180.324952706549742379412290749026850 

- 11339683.06464037751511038174569497508966 e® 

+ 46630964.11789747669532801502683952077042 e® 

- 190369535.4103429881202484621694879170262 + O (e“) 

+ 1.80887954620833474- 12.7814836099524403e 
+ 71.046049240835262 - 334.59648933739741 e® 

+ 1467.5837602507405 - 6165.5621168597119 e® 

+ 25329.619267580422 e® + O (e^) 


(4.12) 


(4.13) 


(4.14) 


Just as in the three-dimensional case m, we have performed various checks on our results. These can 
be divided into two categories: first, we had to make sure that the difference equation, Eq. 12.211 . as well 
as the various parts of its solution, Eqs. itOll . |13|), were in principle correct, and second, that we have 
reached the desired accuracy in the numerical part of our computation. For the first task, it was in general 
enough to ensure that the first few e-orders we obtained for each integral coincided with the existing analytic 
results. Here the main difference to our previous three-dimensional computation laid in the fact that while 
at d = 3 — 2e most of the integrals were either finite or their expansions started with a 1/e term, we now 
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encountered in many cases (often analytically calculable) divergent terms up to 1/e^ order. The analytic 
results relevant to our graphs that we have found in the literature, as well as a few new ones, are collected 
in Sectional We have found agreement in all cases. 

The comparisons with existing analytic results also provide an easy and reliable method to inspect the 
accuracy of the numerical results, since the number of correct digits usually stays roughly constant when 
moving from one e-order to the next. Just as in our previous work |1 Oj . other methods we have employed 
to assess the accuracy question include comparing the results obtained by raising topologically inequivalent 
lines in a single integral to a higher power and analyzing the convergence properties of the factorial series, 
i.e. checking the stability of our results with respect to varying Smax- The results given in the preceding 
section have been observed to be stable to at least the number of digits shown. 

One might be concerned about the rapid growth with increasing e-orders of most of the coefficients. This 
is, as was pointed out in Ref. [7], caused by poles that the integrals (seen as functions of d) develop near 
d = 4, e.g. at d = 112,?,, etc. It is to be expected that factoring out the first few of these nearby poles in 
each case will improve the apparent convergence in e considerably. 

In principle, having a method at hand that is capable of generating coefficients to very high accuracy, 
even to a couple of hundred digits, one could now use the algorithm PSLQ m combined with an educated 
guess of the number content of some of the yet-unknown constant terms, in order to search for analytic 
representations of the numerical results. We have not made any systematic attempts in that direction, since 
the numerical accuracy of our results should be sufficient for all practical purposes. 


5 Laplace transform 


As already mentioned in the above, we have encountered one case where the method of computing the e- 
expansion via a factorial series representation does not work (or, more precisely, does not converge), namely 
for the 3-loop integral of Eq. (Oil . Let us take this specific example as an opportunity to finally display 
a difference equation like Eq. in full detail, and exhibit, following Ref. 0, one method other than 

factorial series for solving it. 

Defining the integral 


M 2 (a;) 




2"~^r(i) 

r(^)r(f) ’ 


(5.1) 


where the dot with the label x means that the corresponding propagator is raised to the x-th power, the 
difference equation Eq. (E3 it satisfies is of second order and reads 


0 = -2{x + l)M 2 {x + 2) + 3{x + 2-d/2)M2{x + l)-{x + 3-d)M2{x) 

r(a: + 5-f) 3-d r{x + 3-d) 1 

r(x + i) r( 5 -f) ’ r(a: + i) r( 2 -d) 

T{x + 2-^) 2 r(a;-b5-M)r(a;-b3-(i) 2 

r(:r + i) r(i-f)^ T{x)r{x + 7-2d) r(i- f) ’ 

with boundary conditions M 2 (a: 1) ^ x~i (cf. eg 112.(ill 1 and 



r(|i)r(i-f)r(f-i)r(|) 

T{d)T{d-2)T{l-d) 


M2(0) 


(5.3) 










We would like to know the master integral M 2 (l), or at least its e-expansion in d = 4 — 2e dimensions. 
Note that only the first two terms of that expansion are known, cf. Eq. below. Formally, it is of 

course possible to solve Eq. (15.21) in terms of factorial series, following the recipe sketched in Section 
However, it turns out that the series does not converge in this (and only this, of all cases treated in this 
paper) case, such that in practice a different method of solving the difference equation is needed. 

One way to tackle Eq. (15.21) could be the iterative method used e.g. in Ref. 52] (cf. Eq. (ll)ff therein): 
expand in e, make Ansatze for the e-coefhcients of M 2 in terms of (sums of multiple) harmonic sums with 
unknown constants, rewrite all e-expansions of the Gamma functions in terms of harmonic sums, then rewrite 
everything in terms of a unique basis, and finally fix the constants by comparing coefficients. Unfortunately, 
there does not seem to exist an algorithm yet that automatizes the choice of Ansatz, hence requiring a fair 
amount of hand-work. For the basic literature on harmonic sums, see the references of Ref. m 

Another way of tackling Eq. 15.21) is to transform it to a differential equation, which should then be solved 
by analytical or numerical methods, or by a combination of both. This is what we will do in the following, 
and this is how we have obtained the numerical values given in Eq. EH). 

Following Ref. 0, we can Laplace transform the difference equation for M 2 , making the Ansatz M 2 {x) = 
/p into a first order differential equation ^o{t)v{t) — = w{t), where 4>o(t) = 3 — d — 

3(1 — |)t — 2t^ and $i(<) = (1 — t){l — 2t). 

The homogeneous equation is solved by vnit) = which however makes 

M 2 (a: 1 ) grow too fast at large x (it would grow like x^~^, in conflict with the large-a: boundary 

condition), such that c// = 0 and hence M 2 ^(x) = 0. 

For solving the inhomogeneous equation, note that the inhomogeneous piece has four terms w{z) = 
which correspond to the last four terms of Eq. (15.21) . written as Tj{x) = dz 
For j = 1,2,3 we therefore have Wj(z) = p^ 3 ^z“^'(l — 2 )““^ where a = (5 — ^,3 — d,2 — |) and 
b = r( 2 -d) ’ ~ r(i-i) )- '^ 4 , we know that it satisfies 


dz z^ ^W 4 {z) = T 4 {x) = 


r(x-b5-^)r(a;-b3-d) 2 

r(a:)r(a:-b 7 - 2d) r(l - |) 


(5.4) 


For T 4 (a:), there obviously is a simple difference equation, x(a:-|-7— 2 d)T 4 (a;-|-l) = (x-|-5—^)(x-|-3—d)T 4 (x), 
from which we get - in complete analogy to Laplace transforming Eq. (15.21) - a differential equation for ^ 4 : 

0 = {d — 3){3d — 10 — Az)w 4 {z) — z{lA — 5d + 4:z{d — 2))w',^{z) + 2z‘^{l — z)w'l{z) . (5.5) 


To write its boundary condition Eq. (Ol) in a (for numerical treatment) more useful form, note that the 
behavior of W 4 {z) at the singular point z = 1 is connected to the large-x limit of T 4 (x). Using Stirling’s 
formula to write = x““^(l -I- _|_ q (x“^) ), we can fix the three constants in the Ansatz 

W 4 {z « 1 ) = ci(l —z)'’’^(l-|-C3(l —z)-|-...), when comparing fg dz z^~^W 4 (z « 1 ) at large x with T 4 (x » 1 ). 
Hence, writing 


, 4(1-z)^-^z^-^W4(z) 

W 4 (z) = Ci(l-zp- - -W 4 (z) = ' ' 


2 sin 


Trd 


d-2 


r(i-|)r(f-i)(d- 2 ) 


■z^-^(l-z)^-^W4(z) , (5.6) 


we get simple boundary conditions ^4(1) = 1^4(1) = ^^(^ 7 ^ — C 3 ) = 2 ^^ ^°r the new function 

W 4 (z), which satisfies the differential equation 


0 = —(d — 4)^W4(z) + z(10 — 3d + 4z(d—3))w',j^(z) + 2z^(z—l)w'J(z). (5.7) 


We now get the non-homogeneous solution v^nit) by varying the constant of the homogeneous solution. 
Due to the linearity of the differential equation, the full solution is simply the sum of four terms, which 
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when plugged back into the definition of the Laplace transform gives a representation for the master M 2 : 


M2{x) = [ zy-^l-2zy-i X 

Jo Jt 

X . (5.8) 


The integral converges (in 4d) for x > 1, so one can use it to compute M 2 ( 2 ) and get M 2 (l) via Eq. (15.211 . 

Unable to solve Eq. (foil for generic d, let us now go to d = 4 — 2e dimensions and start expanding. First, 
we need to solve the differential equation Eq. IS3. Writing W 4 ^{z) = fn{z), the boundary conditions 

translate into /o(l) = 1, /i(l) = -1, /n>i(l) = 0 and /^(l) = 0, /((I) = 0, / 2 ( 1 ) = 2, /),> 2 ( 1 ) = 0. The 
fn{z) satisfy the differential equations 

0 = z(z-l)/"( 2 ;) + (2z-l)/;(z) + (3-4z)/;_i(z) --/„_ 2 ( 2 :) , (5.9) 

Z 

which have to be solved starting with n = 0 (setting fn<o{z) = 0 ). 

One can e.g. solve Eq. 115.911 in terms of multiple integrals. The Ansatz fn{z) = Sn^o ~ ^n,i + fi dag„(a) 
respects the boundary conditions for /n(l) and transforms Eq. (15.911 into a first order differential equation 
for ( 7 „(a), whose boundary conditions gn(l) = 2 Sn ,2 incorporates those for fXO- The homogeneous solution 
is of the form g^io) = aii^a) vanishes due to the boundary condition: c„ = 0. The inhomogeneous 
solution now follows by variation of the constant, such that finally 


fn{z) 

K{a) 


I ~ dn.l + 


da hn{a) 


2(d„,3 - <5n,2)ln(a) 



db (3 — Ab)hn-i{b) 

T 



r ^/i„_2(c) 

h c 


1 — c 


(5.10) 

(5.11) 


The strategy is now clear: /i„(a) —> fn{z) —> W 4 {z) w{z) —> M 2 ( 2 ) ^ ^ 2 ( 1 ). All of these steps can be 

done numerically, and there is a discussion of practical methods in Ref. |3| . 

In practice, we numerically solved for the /„ using Mathematica, changed the order of integrations in 
Eq. (15.811 . dealt with the t-integration (semi-) analytically, and finally performed the z-integration numer¬ 
ically. The singular point at z = 1/2 was treated as a pricipal value integral, and the logarithmically 
divergent regions near z = 1/2 and z = 1 were split off and treated analytically via a series-expansion in z. 

To check the setup, it is possible to start analytically. Solving Eq. (ICTTl . the first couple of orders 
for hn read ho(a) = 0, hi{a) = 0, h 2 {a) = —21n(a) and / 13 (a) = 21n(o) — 31n^(a) -I- 2 Li 2 (l — a), where 
Li„(z) = Yl’k=i W polylogarithm. 

Using Eq. 15.1011 . this then implies fo{z) = 1, /i(z) = —1, f 2 {z) = —In^z — 2 Li 2 (l — z), fsiz) = 
(1-1-5 ln(l — z) — ln(z)) In^(z) -|- 2(1 -|- ln(z))Li 2 (l — z) -|- 10 ln(z)Li 2 (z) — 2 Li 3 (l — z) — 10 Li 3 (z) -|- IOC 3 . 

Knowing now W4(z) = 1 — e -|- e^/2(z) -|- O (e^) and using Eq. 15.611 . we can expand the curly bracket of 
Eq. 15.811 . The two leading terms cancel, such that {..} = ^[(z — l)(7r^ + 21n^(l — z) — 61nzln(l — z)) -I- 
3zln"z + 6Li2(l-z)] + 0(e4) . Now M 2 {x) = dtt^-^ // dz (i_ j(i_2z) (e^) = fo (i-U(i-2U Ci + 

O (e'^) . We obtain M2(2) = 6^36^ -I- O (e"^) , which, using Eq. 15.211 at a: = 0, translates into M2(l) = 
3(4'^^) M2(2) = ^M2(2) = 4^36^ -I- O (e^) , in nice agreement with the first term of Eq. 16.3211 . 


6 Analytic results 

For completeness we list here all existing analytic results applicable to our integrals that we are aware of. 
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Here, we normalize every integral with the appropriate power of the 1-loop tadpole, such that analytic 
results are independent of the integration measure. Also, recall that we have set m = 1. 

We will use the following transcendentals: 

CX5 ^ 

Cn = X] ^ 

k^l 
oo ^ 

o-n = X] 2^^ ^ Lin(l/2) , (6.2) 

k^l 

T 

Lsj(d) = — drln-^”^ |2sin — I , (6.3) 

Jq 2 

and abbreviate the log-sine integrals at their maximum value as Lsj(^) = Ls^ below. 


6.1 1-loop 

There is one 1-loop topology and one coloring by mass. The 1-loop tadpole has an analytic solution in terms 
of Gamma functions. With measure f d‘^p^ J = f d'^P p'^+i = 7r‘^/^r(l — d/2). 

Q = J. (6.4) 


6.2 2-loop 


There is one 2-loop topology and three colorings by mass. One of them reduces to simpler cases, while the 
other two are master integrals. One of the two master integrals has an analytic solution in terms of Gamma 
functions. The other (fully massive) one can be written in terms of the hypergeometric function (see 
Eqs. (4.12) and (4.13) in Ref. JHDj or alternatively in terms of a one-dimensional integral (see Eqs. (21), 
(15) and (16) in Ref. ^^1) which has a simple e-expansion (for 4d in terms of log-sine integrals). 



d—n—2e 


d-2 

2{d-S) 



(6.5) 


r(^)r(i) 

2''-2r(i) 


( 6 . 6 ) 


3(d-2) [ ^ /4-d 5-d 3\ 27rr(5 - d) ] 

4(d-3) ' r(^)r(^)j 


3(d- 2) 
4(d-3) 


a — o , 

l-3~(d-4) 



dr(2sin(r))"^ 


- 3' 


27rr(5 -d) \ 

r(4^)r(^) / 


3(n-2-2e) / , ^ j 4 -n) 3^27rr(5 - n + 2e) | 

4(n-3-2e)j 3 ^ r(4^ + e)r(^ + e) J 


(6.7) 

( 6 . 8 ) 

(6.9) 


The numbers Ls^“^ ~ ~ fo^ dT(2 sin !■)“ In^ ^|2sin^| in the 4d (n = 4) case are the log-sine integrals 
Lsj°^ = Lsj = Lsj (^) of Eq. (16.311 . 
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6.3 3-loop 

There are three 3-loop topologies. 

3-loop, 4 lines: There are four colorings by mass, all of which are master integrals. Two of them have 
an analytic solution in terms of Gamma functions. The third one (called 1,0,1,1,1) in the literature, 

according to the notation introduced in Ref. ini) can be written in terms of a single hypergeometric 
function 3 F 2 (see Eq. (4.33) of Ref. [201 , where also the first seven orders of its 4d e-expansion were given 
in Eq. (4.32))^. The first seven orders of the 4d e-expansion of the fourth (fully massive) master (called 
B^iO, 0,1,1,1,1) in the literature) can be deduced from the function B 4 introduced in Ref. m using the 
reductions Eqs. and given below. Two more orders could be obtained from B 4 as given in 

Ref. |24| . but we refrain from reproducing them here. 



(ID 




r3(^) 

2^-3r(^)r(^)r(f) 


( 6 . 10 ) 

( 6 . 11 ) 


d=4-2e 


~ ~ ^v^Ls 2 ^ -I- - ^D3(7r^ -f 9(3-21n3)Ls2 -flSLss)^ 

/ TQOT 1 -7T- 

-f -f -D3(97rD21n3 - 3) -f 64Ls4(-) - 9(67 - 54In3 -f 18 In^ 3)Ls2 

y 64 8 3 

-H 162(2 In 3 - 3)Ls3 - 2 I 6 LS 4 -f 1847rC3) 


/48601 

I 128 


( 3 23 69 

— 7r^(67 -I- 181n3(ln3 — 3)) — -|- 7r(3 — 21n3)^3 -I- 8l7rLs( 

J-0 oO Z 

-f -^(-457 -f 6 In 3(67 -f 3 In 3(2 In 3 - 9)))Ls2 - — (67 -f Stt^ -h 18 In 3(ln 3 - 3))Ls3 
16 8 

+ -^(2 In 3 — 3)Ls 4-+ 12(3 — 2 In 3)Ls4( — ) -I- 49Ls5( —)-^Ls(^ ^ e® 

+ 0(eD (6.12) 


d=4-2e _ 5 1 2 103 3 ^ 

3 2 12 24 


(163- 128C3)e'‘ 


( 


9055 
/63517 


——-h In^ 2(7r^ — In^ 2) — I68C3 — 25604^ 

4o o / 


-\ -In"^ 2(4 In2 — 15)- tt^ In^ 2(4 In2 — 9)-7r'^(41n2 — 3) 

V 96 5 ^ 3 ^ 15 ^ ’ 

+ I 24 OC 5 - 115204 - 153605 ) e® + O (e^ 


(6.13) 


Here, Ls' = — Jq^ drT^ ^ In^ |2 sin || are special values of the generalized log-sine function ^Hj . 

^In some sense, the representation in terms of special types of hypergeometric functions can be called an all-order analytic 
e-expansion, namely when their expansion can be written in terms of rapidly converging (multiple inverse binomial) sums, for 
which efficient algorithms exist m The 3 -loop integrals E3, D5 and D4 below belong to this class as well. 
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3-loop, 5 lines: There are eleven colorings by mass. Eight of them reduce, while the remaining three 
are master integrals. One of the masters has an analytic solution in terms of Gamma functions. The second 
one (called E 3 in the literature) can be written in terms of the hypergeometric function 2 F 1 , cf. Eq. (4.24) 
of Ref. 120]. Its first six terms of the 4d e-expansion (we will only reproduce the hrst five of them below) 
are given in Eqs. (4.16),(4.18) of Ref. jOO]- The first five terms of the 4d e-expansion of the third (fully 
massive) master integral (called in the literature) can be deduced from Ref. m using the reduction 
given in Eq. (16.2511 below. One more term has recently been given in Eq. (3.28) of Ref. |26| . but we refrain 
from listing it here. 




J3 


J3 




J3 


QD - 

6(d-3){‘“ ®*(D./ ^>© 0 } 

(6.14) 



(6.15) 

QJ) = 


(6.16) 


-5ibi|»-aD-^(o)‘) 

(6.17) 


®©)-yeo-i®(o)’ 

(6.18) 


3d - 8 


4(2d- 7) 


d-2 
3d - 8 
“ d-4 


r2(<i^) 


=4-2e 2 

3 


siiR(f) sin(^) r3(^)r2(d - 2 )r(f) 

-b -b (5 -b y - 3 V 3 Ls 2 ^ 


(6.19) 

( 6 . 20 ) 
( 6 . 21 ) 

( 6 . 22 ) 


(6.23) 


+ (y + y + (S 2 + ) '' 

+ —I —Q -y*^^ ^^^33(10 -b (ln3 — 4) ln3))Ls2 

-b3(ln3 — 2 )Ls3 — 2Ls4 — ^^® 4 (-^) + — 2) -b -b O (e^) 

1 -b -b - 6y3Ls2^ -b - 6C3 + ^3 + 6(ln3 - 2)Ls2 - 6LS3 

-b (76 - ^ -b I8LS2 - 127rLs3 -b I8LS4 + (^-y + 4\/37r -b 26 In 3^ Ca 
-bVs ^y(ln3-2)-3(10-41n3-bln2 3)Ls2-b6(ln3-2)Ls3-4Ls4^^ e^ + O (e^) (6.24) 
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3-loop, 6 lines: There are ten colorings by mass. The first two terms of their 4d e-expansion are given in 
Ref. |25|. Five of them reduce, while the remaining five are masters. The third term of the 4d e-expansion 
of one of the masters (called D4 in the literature) is given in Eq. (4.10) of Ref. m^- The remaining four 
masters (called D^, £>3 and Dq, respectively) are read from Ref. m 



© ^ 



J3 



J3 




J3 


2(d-3) /T7\ 3d-8 2(d-2)^ ^ 

3{d-4)K\U l2{d-A)\SJJ 3(d-4)Ow 


(d-2)^ 


6(d- 4)(d- 3) 
(3d-10)(3d-8) 


16(d-4)2 

3(3d- 10)(3d- 8) 
' 16(d-4)(2(i-7) 

(3(i- 10)(3d- 8) 


o 

(ID 


(d-4)2 
2(d- 3) 
d-4 


4(d-4) 

2d-7 viD.. 
{d-2f 




8(d-4)((i-3) 
d-4 


(d-2)2(5d-18) 

8(d-4)2((i-3) 

Q ’ 


o 


4(2d-7) 
2{M- 10) (3d-8) 


d-2 

d-4 


eo 


(6.25) 

(6.26) 

(6.27) 

(6.28) 
(6.29) 


-2C3e" 




21 

—5-X5 + 


(1080 
161 


2 ^ 


7rLs4 (—) 
54 ^^3' 


367 3, „ ^ 2615 2. 2047, 

-7r3Ls2 — 77rLs4 — 2 C 3 H--Cs 

216 ^ 432 ^ 216 ^ 


0(eD 


(6.30) 


- 2 C 3 D -f ^ + 9Ls2^ + 0 (e^) (6.31) 

- 2 C 3 e^ + + I iD 2(^2 _ 2) _ ^3 ^ 33) 

- 2 ( 36 ^ + + 0 (e^) (6.33) 

-2(36^ + + \ + ®^®2 - 1604^ + O (eD (6.34) 


Here, xs = 7^^;^ Dj=i } ~ 0.0678269619272... is a special case of a binomial sum ^H] . 


6.4 4-loop 

There are ten topologies. 

4-loop QED-type cases, 5 lines: There is one topology, BB. 

There are two QED-type colorings of BB. Both of them are masters. One is known analytically in terms of 
Gamma functions, while the other one is new. Interestingly, the analytic value of the last term in Eq. 

“^Note that there is a typo in Eq. (4.10) of Ref. EOl- The second-last term should read — ^-7rLs4(-^), see also Ref. 1261 . 
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was obtained by a physics computation in which this master integral contributed m 



= 3((i-2)4'^-3 


r(5 - 2d)r(^)r(^)r2(|) 

r(ii^)r3(i^ 


(6.35) 


d=4-2e ,1 17 2 


1 3 37207 4 

~ 1296 ^ 


^^ + ^C3)e^ + 0(e«) (6.36) 


4-loop QED-type cases, 6 lines: There are two topologies, T and G. 

There are four QED-type colorings of T. All of them are masters. One of them is known analytically, 
while the first six orders of the 4d e-expansion of two others were given in Eq. (16) of Ref. [7] and Eq. (18) 
of Ref. respectively. The fourth one is new. 






J4 



J4 


8d-3r3(i)r(6-2d)r3(f) 

sin(3|d)r(ii^)r2(4^)r2(d - 2 ) 


(6.37) 


3 7 9 2 

2+ 2"+2' 


-1254- 




-208+ —+ 109C3 )e^ 


5477r4 

60 


+ 32 In^ 2(ln2 2 - tt^) + 76804 + 855C3 + ISOCs e' 


+ O (e®) 


(6.38) 


3 3*^ 


+ - 


+ O (e^) 


1928 3267r4 64. 2,,„„2o 2^ , , 1192C3 


H—— In 2(ln 2 — tt ) + 51204 + 

45 3 3 


+ 96C5 e' 


(6.39) 


1 1 


+ -e + 0 • e^ + -8 + —C3 e 


13 


241 57r4 


+ 4C3 e 


+ - 


669 TT^ 


+ 36C3+“c 5) e=+0(<!«) 


(6.40) 


There are five QED-type colorings of G. All of them reduce. 



2 d -5 /7TV^ 
~6(d-3) viLiy 



(6.41) 

(6.42) 

(6.43) 


^Note that in Ref. m the last term of Eq. 16.391 involves a numerical coefficient A^io ~ 5.3111546, which we have determined 
to be Nio = ^ ^ ~ ~ 4a4, using our high-precision result Eq. 14.101 and PSLQ 1161 . 
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(2d-5)(3d-8) 

“ 6(d-3)(d-4) viiZy 

2d-5 

2{d-3) vJJiy ~ 2(fi-3) 


S)0 


(6.44) 

(6.45) 


4-loop QED-type cases, 7 lines: There are three topologies, VB, N and U. 

There are seven QED-type colorings of VB. Five of them reduce. There are two master integrals. Both 
are new. 





(© 


(2d-5)(3(i-10)(3d-8) /TTTN , 2{d-3f 

3((i-4)2(3d-11) vAliy ^ 

(2d - 5)(3d- 8) 

3(d-4)2 'oiiy 


(d-4)(3d- 11) 't; 
4(d-3)2 


3(d-4)(3d- 10) 


(d-2)(3d-8) ^ 

8{2d - 7)(3d - 10) Lv 

(2d - 5) (3d - 8 ) (rf - 3) ,■ 

3(d-4)(d-3) vjJZy ~ d- 4 
(2d-5)(3d- 8 ) 2(d-3) 


16(d-4)(d-3) 

(d-2)3 


3(d-4) 


32(d-4)(d-3)2 
2 


o 


3d-10 


(d-3)(^ - 


3d - 8 
2(d-3) 


(6.46) 


(6.47) 




(d-2)(3d-8) ^ 

8(d-4)(d-3) vAZy W 


^(D 


d -2 


( i ))0 


(6.49) 

(6.50) 


d=4-2e 1 5 , 


11 

3" 

.4 


44 y 2.3 
3 60 3^^ ^ ^ 


332 TT^ 31 


+ O (e^) 


(6.51) 



J4 


rf=4-2e 



(6.52) 


There are five QED-type colorings of N. All of them reduce. 



(2d-5)(3d- 8 ) ^7TT^ (d-2)(3d- 8 ) {d-2)^ 

6(d-3)2 ZZZy 4(d-3)2 W " 8(d-3)3 

(2d-5)(3d-8) ^7rr\ (d-2)(3d-8) ZTS"'-. 

" 6(d-4)(d-3) VZLiy ^ 8(d-3)(2d-7) vAZ.-'' W 
(2d-5)(3d-8) ZTlZy 

6(d-4)(3d- 11 ) vAlAy 

(2d-5)(3d-8) (d-2)(3d-8) ^ 

i6(d-4) vAZy 4(2d-7) vAZ.yw 


(O )*<•=■“' 

(6.54) 

(6.55) 


3(d- 2)3 
32(d-4)(d-3) 



(6.56) 
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(6.57) 



(2d-5)(3rf-8) ^^ {d-2){3d-8) 

16(d-3) 8(d-4)(d-3)ViZyW 

_ (d-2)(3rf-8) ^ , (d-2)3(3d-4) 

8{d-3) vi,/..-''VJ 32{d-4){d-3)^ J 


There are four QED-type colorings of U. All of them reduce. 




d-3 
d — 5 
d-3 
'2(d-4) 

d-3 
'3d-11 


id-2){3d-8) ^ _ (d-2)3 

8(d-3)2 8(d-3)3 

Vs 3(d-2)(3d-8) ^7Y^ 3(d-2)3 

8(d-5)(d-4) VW W 4(d-5)(d-4)(d-3) 
, 7 ^ _ (d-2)(3d-8) ^ 

'O' 8 (d - 4)(2d - 7) W...-' W 



(6.58) 

(0)*'““' 

(6.60) 

(6.61) 


4-loop QED-type cases, 8 lines: There are two topologies, VV and W. 
There are seven QED-type colorings of VV. All of them reduce. 



(2d-7)(2d-5)(3d-10)(3d-8) (d-3)2(2d-7) 

6(d-4)2(3d- 13)(3d- 11) OiV ~ (d-4)(3d- 13)(3d- 11) 'V 
(2d-7)(2d-5)(3d-10)(3d- 8 ) (d-3)(2d-7) y:’:;-, 

18(d-4)2(d-3)(3d-ll) Viliy ~ 3(d-4)(3d- 11) 'O' 

(d - 2) (3d - 10) (3d - 8 ) 

~ 48(d-4)(d-3)(2d-7) O'..--' O' 

(2d-7)(2d-5)(3d-10)(3d-8) (d-3)(2d-7) .-'/’Vs 

6(d-4)2(d-3)(3d- 11) Oiy ^ 4(d-4)2 'O' 

(d-2)(3d-8)(5d2-35d-f 61) O) 

^ 16(d-4)2(d-3)(2d-7) O'..--' O' 



{2d - 7)(2d - 5)(3d - 10)(3d - 8 ) OTO 

^ - 

6(d-4)2(d-3)(3d-11) OZy 


4(d- 3)2(2d- 7) 
3(d-4)(3d- ll)(3d- 10) 


(d - 2)(3d - 8)(95d3 - 989d^ + 3428d - 3956) (d - 2)^ 

32(d-4)(d-3)(2d-7)(3d- ll)(3d- 10 ) O'..--' O/ 16(d- 4)(d- 3)2 


(6.62) 


(6.63) 


(6.64) 


(6.65) 

( 0 )‘ 


AZZ 

(2d-5)(3d-ll)(3d-8) 

07 " 

32(d-4)(2d-9) OLZ 


(d - 2)(3d - 8)(12d2 - lOld-f 204) 
32(d-4)(2d-9)(2d-7) 


id-3r 

4(d-4)(2d-9) 


GDO 


(3d-10) 

4(2d-9) Oiy 


(d- 2)(3d- 8 ) 

~ 16(d-4)(2d-9) 
(2d-5)(3d-8) 
12(d-3)(2d-7) OZy 


ODO 


3(d- 2)3(3d- 7) 


64(d-4)(d-3)(2d-9) 
(2d - 5)(3d- 8 ) 


64(d- 3)(2d- 7) 


o 


2(d-3)2 


3(2d- 7)(3d- 10) 'O'' 


3d-10 /'V-Oi (d-2)3(73d2-512d-f 896) 

2 ( 2 d- 7) '■•-..'O' ~ 128(d-4)2(d-3)2(2d-7) V O^ j 
(d-2)(3d-8)(19d2-128d-f 216) OS”''-- O) 

~ (16(d - 4)(d - 3)(2d - 7)(3d - 10 ) Ol'..--' O/ 


( 6 . 66 ) 


(6.67) 
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(d- 2)(3d- 10)(3d- 8) 


32(d-4)2(d-3) 

(2d - 7)(2d - 5)(3d - 10)(3d - 8) 


(®o 


64(d-4)2(d-3) 
(d-2)(3d-10)(3d-8) 
8(d-4)(2d-7) 
(d-2)3(18d3- 129d2- 


+ 


(d-2)3(73d2-512d + 896) 
128(d-4)2(d-3)2(2d-7) 
(d- 3)(2d- 7) 
3(d-5)(d-4) 


o 


2)0 


(d - 2)(3d - 8){lld2 - 77d + 134) 


32(d-5)(d-4)2(d-3) 


(OO 


245d- 58 


128(d-5)(d-4)2(d-3)2 


o 


( 6 . 68 ) 


There are five QED-type colorings of W. Four of them reduce. There is one master integral, which is new. 


0 


(2d - 5)(3d - 8)(27d3 - 2S2,(P + 990d - 1156) (d - 3)(3d - 10) 


+ 


12(d-4)3(2d- 7)(3d- 11) 
2 (d- 3)25d- 18) 

3(d - 4)2(2d - 7)(3d - ll)(3d - 10) 

(d- 2)(3d- 8)(5d- 18) 
16(d-4)(3d- ll)(3d- 10) 


OJO/ 2(d-4)(2d-7) 
(2d-5) (3d-8) 


2)0 


64(d-4)(2d-7) 
7(d- 2)3 


(6.69) 


CB 


2(2d - 7)(2d - 5)(3d - 10)(3d (^j)^ 


+ 


9(d-4)3(3d- 11) 
2 (d-3)3(2d-7) 


128(d-4)(d-3)(2d-7) 
8 (d- 3)3(2d- 7) 


o 


9(d-4)2(3d-ll)(3d-10) '00' 
(d - 2) (3d - 8) (7d2 - 48d + 82) 


3(d-4)2(3d- 11) ■'O 
(2d-5)(3d-8) OOO 

6(d-4)(2d-7) 'BLiy 

4(d- 3)3 

~ 3(d-4)(2d-7)(3d- 10) 

(d- 2)(3d- 8) 


24(d-4)(2d-7)(3d-ll)(3d-10) 
(2d-5)(3d-8) . d-3 


2)0 


(6.70) 


32(d- 3)(2d- 7) 

(d-3)(3d-10) , 
(d-4)(2d- 7) ' 
(d- 2)(3d- 8) 


2(d-4) 'y 

3d - 10 


+ 


8(2d- 7)(3d- 10) 
(d- 2)3(9d- 28) 


0 


64(d-4)(d-3)2(2d-7) 
(2d - 5)(3d- 8) 


0)0 

o 


16(d-4)(d-3) 


2(d-4) CS) 

(IDO 


(6.71) 


+ 


+ 


32(d-4)(d-3) 
3d - 10 


+ 


2(d- 3)2(2d- 7) 


3(d-4)2(3d- 11) 
(d-2)(3d-8)(d2-4d + 2) 
2(d-4) ' 8(d-4)2(d-3)(3d- 11) 

(d-2)3(29d2 - 177d + 268) ^ ^ 
64(d-4)2(d-3)2(3d- 11) 


d-3 


(O 


2(d-4) 

(OO 


o 


d=4-2e 


J4 


5C5e^ + 0(e") 


(6.72) 


(6.73) 


4-loop QED-type cases, 9 lines: There are two topologies, H and X. 
There are five QED-type colorings of H. All of them reduce. 


a 


D ^ - 


3(2d - 5)(3d - ll)(3d - 8)(54 - 29d -f 4cP) 

256(d - 5)(d - 4)2(d - 3)(2d - 9) 

3(d-3)(d2 - lld-H27) ,-vV\ 9(d - 4)(3d - 10) 


8(d-5)(d-4)(2d-9) 


8(d-5)(2d-9) 


4(d-3)y2d-7) 

~ 3(d- 5)(3d- 13)(3d- 11) 

(O 
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id-2){3d-8){13d^-77d^+ 9d + 351) {"7^^ 

32(d - 5)(d - 3)(2d - 9)(3d - 13)(3d - 11) vlZy 
3(d - 2)(3d - 8)(13d2 - 92d + 162) 

64{d - 5)(d - 4)(2(i - 9){2d - 7) vW..--' 


(6.74) 


(d - 2)3(2452^5 _ 43031^4 + 329345^3 - 1198763^2 + 2170827d - 1564110) 
512(d - 5)(d - 4)2(d - 3)2(2d - 9)(3d - 13)(3d - 11) 


o 


6 (d- 3)2(2d- 7) 

(3d - 14)(3d- 13)(3d- 11) 


3(2d - 7)(2d - 5)(3d - 10)(3d - 8) 

“ 2(d - 4)(3d - 14)(3d - 13)(3d - 11) 

32(d- 3)3(2d- 7) 

~ 9(3d - 14)(3d - 13)(3d - ll)(3d - 10) 

(d - 2)(3d - 8)(139d3 - 1495d2 + 5344d - 6348) 77S''\ 

96(2d - 7)(3d - 14)(3d - 13)(3d - ll)(3d - 10) viJ..-'' 

7(2d-7)(2d-5)(3d-10)(3d-8) ^7)^ , 4(d- 3)2(2d- 7) 

9(d-4)2(3d-13)(3d-ll) vjJZy 
32(d- 3)3(2d- 7) 


(6.75) 


3(d-4)(3d- 13)(3d- 11) 
d-4 


9(d - 4)(3d - 13)(3d - ll)(3d - 10) 2(3d - 13) 

(d - 2)(3d - 8)(409d3 - 4285d2 + 14944d - 17348) 


96(d - 4)(2d - 7)(3d - 13)(3d - ll)(3d - 10) 

4 


070 


(d-2)3 


o 


16(d-4)(d-3)(3d-13) 

(2d - 5)(3d- 8)(3d3 - 19d^ 

12(d - 5)(d - 4)2(2d - 7)(3d - 11) 
(2d - 5)(3d - ll)(3d - 8)(2d3 - 14d' 


(6.76) 


19d+38) (^7^ 
23d + 6 ) 


128(d - 5)(d - 4)3(d - 3)(2d - 7) 

2(d-3)2(2d-7) /7^ 2(d-3)3(17d2-125d+ 230) 


3(d- 5)(d-4)(3d- 11) 

(d-3)(5d2 - 37d + 69) ,-■/ 

4(d-5)(d-4)2 ’ 2(d-5)(2d-7) 

(d - 2)(3d - 8)(10d3 - 87d‘^ + 235d - 186) 


3(d - 5)(d - 4)(2d - 7)(3d - ll)(3d - 10) 
3(d-3)(3d-10) 3(3d-10) 


4(d- 5) 


(© 


GDO 


GDO 


32(d-5)(d-4)2(d-3)(3d- 11) 

(d - 2)(3d - 8)(38d3 - 568d3 + 3174d2 - 7863d + 7290) 

16(d - 5)(d - 4)2(2d - 7)(3d - ll)(3d - 10) 

(d - 2)3(418d5 - 7346d4 + 51389d3 - 178846d2 + 309603d - 213234) 
256(d - 5)(d - 4)3(d - 3)2(2d - 7)(3d - 11) 

(2d - 5) (3d - 8)(69d3 - 725d^ - 2543d - 2978) (^7^ 


o 


(6.77) 


12(d - 4)2(2d - 9)(2d - 7)(3d - 11) 
3(2d- 5)(3d- 8 ) 


128(d- 3)(2d- 9)(2d- 7) 
(d-3)(d2-lld + 27) /;■ 
8(d-4)2(2d-9) 

3(3d- 10) 


2(d-3)3(d-2)(13d-47) 


9(d - 4)(2d - 9)(2d - 7)(3d - ll)(3d - 10) '77' 
3(d-3)(3d-10) .•■' 


8(2d- 9) 


(® 


+ 


2(2d-9)(2d-7) '-..V 
(d- 2)(3d- 8)(5d- 18) 


64(d-4)2(d-3)(2d-9) 


GDO 
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{d - 2){3d - 8)(1139(i4 - 16453^3 + 89068^2 - 214178 + 193044) 
192{d - 4)2(2d - 9){2d - 7)(3d - ll)(3(i - 10) 

{d - 2)3(48^3 - 445^2 + 1352d - 1344) ^ ^ ^ 


02)0 


O 


256(fi - 4)2(d - 3)2(2(i - 9)(2(i - 7) 

There are two QED-type colorings of X. One reduces. The other one is a master integral. It is new. 

(2d - 5) (3d - 8) (2109d'‘ - 31288d3 + 173302d2 - 425005d + 389562) 


(6.78) 


Cfi; 


(H) 


J4 


+ 


+ 


+ 


+ 


36(d - 4)(2d - 9)2(2d - 7)(3d - 13)(3d - 11) 
(2d - 5) (3d - 8)(24d3 - 268d'^ + 1003d - 1258) 
256(d-4)2(2d-9)2(2d-7) 

14(d-3)2(2d-7) //V-. 


3(2d-9)(3d-13)(3d-11) 

2(d - 3)3(295d3 - 3332d2 + 12431d - 15334) 
9(2d - 9)2(2d - 7)(3d - 13)(3d - ll)(3d - 10) 
3(d-4)(d-3)(3d-10) 


2(2d-9)2(2d-7) 


+ 


2(2d- 9)(3d- 13) 


+ 


(d - 2) (3d - 8)(599d3 - 9067d3 + 51340d2 - 12886d + 121044) 
48(2d- 9)2(2d- 7)(3d- 13)(3d- ll)(3d- 10) 

(d - 2)3(392d4 - 6204d3 + 36843d2 - 97323d + 96502) 


ODO 


512(d - 4)2(d - 3)(2d - 9)2(2d - 7)(3d - 13) 


o 


d=4-2€ 


Xoe^ + O (e^) 


(6.79) 


(6.80) 


All the above formulas agree with our numerical results of Section^ 


7 Conclusions 

We have employed the general method of numerically solving single-scale integrals in terms of their e- 
expansion around d = 4—2e via difference equations, to high precision and to high e-orders. We have covered 
the set of all vacuum master integrals up to three loops, as well as ‘QED-type’ vacuum master integrals at 
4-loop order. These integrals play a role in state-of-the-art perturbative calculations for precision tests of 
the standard model. 

The main vehicle of solving the difference equations treated in this work was a formal representation in 
terms of factorial series, which could then be evaluated numerically in a truncated form. 

In cases where the factorial series representation does not converge, a more general (and hence more 
complicated) method can be used, which transforms the problem into differential equations. We have 
encountered only one such case, and have shown in detail how it can be represented in terms of multiple 
integrals, which we then solved numerically. 

Furthermore, we have made an attempt to collect all presently known analytic results for the class of 
vacuum master integrals that we have treated here, up to the 4-loop level. This is meant as a concise 
reference for practitioners in the field. 
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A Numerical results for analytically known master integrals 


As a complement to Section 0] we here list the first few terms of the Laurent expansions in e = (4 — (i)/2 
of those single-mass-scale vacuum master integrals up to four loops that are known analytically (see the 
explicit d-dimensional expressions of Section EJ. 

Notation and integral measure are as in Section 01 which in particular determines the 1-loop tadpole to 


be J = 


r(2+^ 


I- 


p2 + l 




E OO 

n=C 


^2n-l 





- 1.0000000000000000000000000000000000000000000000000 £-1 

- 1.0000000000000000000000000000000000000000000000000 e 

- 1.0000000000000000000000000000000000000000000000000 e^ + O (e®) (A.l) 

- 0.5000000000000000000000000000000000000000000000000 £-2 

- 0.5000000000000000000000000000000000000000000000000 £-1 
-3.6449340668482264364724151666460251892189499012068 

- 3.4428771636886321510726770051345751984539636088663 £ 

- 17.748133915933433322311939139507906328597192596121£2 

- 16.366439374126401323669287645924253086404587829687 + O (e"^) (A.2) 

- 0.083333333333333333333333333333333333333333333333333 £“2 

- 0.37500000000000000000000000000000000000000000000000 £-1 

- 2.4683003667574465515695409166563459279428082839367 

- 8.5848042311088475775631523236940150167718153674315 £ 

- 38.120827450450135424466436253406610052456582985006 + O (e^) (A.3) 

-k 0.33333333333333333333333333333333333333333333333333 £"^ 

-k 0.16666666666666666666666666666666666666666666666667 £"2 
-k 0.58333333333333333333333333333333333333333333333333 £"^ 

-k 0.41381840842558476106596843069719997537329677957466 

- 24.905969600320865917659060143145414845610363033237 £ 

- 12.059724940640299353325034075589267393005352211165 + O (e^) (A.4) 

- 0.33333333333333333333333333333333333333333333333333 £"^ 

- 0.66666666666666666666666666666666666666666666666667 £"2 

- 5.9565348003631195396114969999587170451045664690803 e-^ 

- 10.976993729846780032023343117902167435855817881707 

- 67.587197404302297868575437012376235190940093056288 £ 
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- 120.04368176952781664333978568556636327245719008030 + 0 (e^) (A.5) 

+ 0.08333333333333333333333333333333333333333 
+ 0.2361111111111111111111111111111111111111e"^ 

+ 0.4189814814814814814814814814814814814815 e"^ 

+ 0.5870437170675600697079437393391752840153 

- 38.15649063807203021274518890590965524153 e + O (e^) (A.6) 

+ 0.1666666666666666666666666666666666666667 
+ 0.1666666666666666666666666666666666666667 

- 0.1666666666666666666666666666666666666667 
+ 0.1130270533440971408650083562391832748449 

- 78.06045784302350379977045250776444270422 

- 77.68037792722049725359022296531379645497 e + O (e^) (A.7) 


In fact, the results shown here have been obtained via numerically evaluating truncated factorial series 
along the lines of Sections [3 and |3 but they of course coincide perfectly with the analytical results of 
Section 13 (note the different normalization there). 
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